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Abstract 

The entanglement entropy in three-dimensional conformal field theories (CFTs) receives 
a logarithmic contribution characterized by a regulator-independent function a(9 ) when the 
entangling surface contains a sharp corner with opening angle 9. In the limit of a smooth 
surface (0 —>• 7r), this corner contribution vanishes as a{9) = a (9 — tt) 2 . In arXiv: 1505.04804, 
we provided evidence for the conjecture that for any d = 3 CFT, this corner coefficient a is 
determined by C T , the coefficient appearing in the two-point function of the stress tensor. Here, 
we argue that this is an instance of a much more general relation connecting the analogous 
corner coefficient a n appearing in the nth Renyi entropy and the scaling dimension h n of the 
corresponding twist operator. In particular, we find the simple relation h n /a n = (n — l)7r. We 
show how it reduces to our previous result as n —> 1, and explicitly check its validity for free 
scalars and fermions. With this new relation, we show that as n —> 0, a n yields the coefficient 
of the thermal entropy, c s . We also reveal a surprising duality relating the corner coefficients 
of the scalar and the fermion. Further, we use our result to predict a n for holographic CFTs 
dual to four-dimensional Einstein gravity. Our findings generalize to other dimensions, and 
we emphasize the connection to the interval Renyi entropies of d = 2 CFTs. 
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1 Introduction & main results 

Understanding the structure of quantum entanglement in complex systems has become an active 
area of study in a variety of areas of physics, including condensed matter physics, e.g., [1, 2]; 
quantum field theory, e.g., [3, 4]; and quantum gravity, e.g., [5, 6]. For these investigations, the 
entanglement entropy S EE and Renyi entropies S n [7] have proven to be two particularly useful 
measures of the relevant degrees of freedom. In the context of quantum field theory, these are 
defined for a spatial region V with: 

S n (y) = w—— log Tr py , S EE {V) = lim S n (V) = - Tr (p v logp v ) , (1.1) 

where py is the reduced density matrix computed by integrating out the degrees of freedom in the 
complementary region V. 
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Figure 1: a) Region V whose boundary contains a sharp corner with opening angle 9. b) The 
contribution to the Renyi entropy S n from a corner in the almost smooth limit yields a great deal 
of insight into the degrees of freedom of the CFT via the coefficient a n . 


In the present paper, we will focus on the Renyi entropy in three-dimensional conformal field 
theories (CFTs), which takes the form 

Sn = bJ- - a n {9) log {i/8) + c n + 0(8/i) , (1.2) 

where 5 is a UV cut-off and i is a length scale characteristic of the size of the region V — see 
Fig. 1(a). The result is dominated by the first contribution, the celebrated ‘area law’ term, but the 
corresponding coefficient B n depends on the details of the UV regulator. The subleading logarithmic 
contribution appears when the entangling surface (i.e., the boundary of V) contains a sharp corner 
of opening angle 9, as in Fig. 1(a). The corresponding corner function a n {9) is regulator independent 
and hence it is a useful quantity to characterize the underlying CFT. For instance, several groups 
have numerically studied the corner function using lattice Hamiltonians [8, 9], and obtained results 
independent of the lattice details which probe the low energy degrees of freedom. 

As the corner function a n {9) is central to our investigation, let us summarize a few of its key 
properties: If eq. (1.2) is evaluated for the vacuum state, reflection positivity [10, 11] constrains 
a n {9) to be a positive convex function of 9, i.e., 

a n {9 ) > 0, 8 e a n (9) < 0, d 2 0 a n (9) > 0, (1.3) 

in the range 0 < 9 < tt. In fact, reflection positivity gives rise to an infinite tower of nonlinear 
higher-derivative constraints as well [10]. 1 Further, this function satisfies a n {9) = a n (2 ti — 9) if 
eq. (1.2) is evaluated for a pure state, e.g., in the vacuum state of the CFT. The form of a n {9 ) is 
also constrained on general grounds in the limits where the corner becomes very sharp (9 —y 0) and 
where it becomes almost smooth ( 9 —y 7r, Fig. 1(b)) [4]: 

a n (9 -> 0) = y , a n (9 ->• vr) = a n (tt - 9) 2 . (1.4) 

This behaviour is schematically illustrated in Fig. 2(a). Hence these limits define two regulator- 
independent coefficients, K n and a n , which are representative of the CFT. 

In studying the corner contribution in the entanglement entropy [12, 13], we recently conjectured 
that the smooth-corner coefficient o\ has a universal form in general three-dimensional CFTs, 

ai = 24 ’ (1-5) 

1 We thank Horacio Casini for explaining these points. We discuss the nonlinear constraints further in section 6. 


2 








where C T is the central charge appearing in the two-point function of the stress tensor — see 
eq. (1.12). We have verified that this relation holds for a free conformally coupled scalar and a free 
massless fermion, as well as for an eight-parameter family of strongly coupled holographic CFTs 
[12, 13]. A more general holographic proof appears in [14, 15]. Our primary result here is the 
generalization of eq. (1.5) to arbitrary values of the Renyi index n > 0: 




1 h n 
7T n — 1 ’ 


( 1 . 6 ) 


where h n is the scaling dimension of the twist operator r n appearing in calculations of S n (V). It is 
defined in the n-fold replicated theory as the surface operator on the boundary of V which acts to 
permute the n copies of the original QFT — see section 2 for more details. 

A first test of eq. (1.6) is to verify that we recover eq. (1.5) from this new relation in the limit 
n —y 1. To do so, we make two observations: First, the twist operator becomes trivial at n — 1 and 
hence the scaling dimension h n vanishes in this limit. Second, at n — 1 the first derivative of h n 
with respect to n is proportional to the central charge C T in any d -dimensional CFT [16, 17]. In 
d = 3, the precise relation is 


9,ih n \ n=1 



(1.7) 


Therefore in the vicinity of n = 1, the scaling dimension can be expanded as 


( 1 . 8 ) 


Now it is straightforward to see that substituting this expression into eq. (1.6) and taking the limit 
n —>■ 1 precisely reproduces the original relation (1.5). 

Further, we have been able to verify eq. (1.6) for a free conformally coupled scalar field and for 
a free massless fermion. In particular, as we discuss in the following, we can evaluate the corner 
coefficient a n using the results of [4, 18] while the scaling dimension h n can be determined with 
the results of [16, 17]. We demonstrate that these two independent calculations yield complete 
agreement with eq. (1.6) for any integer values of n > 1. Since the methods involved in these two 
calculations are so completely disparate, we find this agreement to be very strong evidence for the 
new conjecture. 

Before moving to detailed discussions, let us draw an interesting parallel with two-dimensional 
CFTs where twist operators are well understood. In the case of d = 2, the twist operator is a local 
primary operator with scaling dimension [3, 19]: 


K ] = — 
n 12 



(1.9) 


where c is the Virasoro central charge of the theory. The Renyi entropy of a single interval is 
calculated by first evaluating the correlator of two twist operators inserted at each of the endpoints 
and the final result can be written as 

/? (2) 

Sn ~ 2 = 2 °n !°g( V 5 ) + ' • • 1 where °n = ^~T • (1-10) 

n — 1 
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Figure 2: Schematics for the corner entanglement function in general CFTs in d = 3 spacetime 
dimensions, a) Renyi corner contribution a n {9) versus the Renyi index n and opening angle 6 , with 
the asymptotics for 6* —0,7r and n —y 0. b) Smooth-corner coefficient cr n versus n. The divergence 
at small n is determined by the flat-space thermal entropy density: s flat = c s T 2 . 


In this expression, t is the length of the interval and S, the UV cut-off. The factor of two here comes 
from having two endpoints [19]. Hence, there is a striking similarity between the coefficient of the 
logarithmic contribution in d — 2 and in d — 3 — in the limit 6 —y n for the latter. The parallel 
extends to the n —> 1 limit, where one recovers the well-known result S^ 2 = (c/3) log (£/S) by 
applying the two-dimensional analog 2 of eq. (1.7): d n h n \ n= \ = c/6. While it calls for some deeper 
physical insight, this interesting connection serves as a stepping stone from the well-known results 
for d = 2 to new considerations of Renyi entropies in higher dimensions — see further discussion in 
sections 4 and 7. 


1.1 Main results 

Eq. (1.6) is our main result. In particular, we conjecture this to be a relation valid for general CFTs 
in d — 3. We thus predict that the corner coefficient of the Renyi entropies in the smooth limit 
<j n is proportional to the scaling dimension h n of the twist operator. Most of the remainder of the 
paper is devoted to supporting this new conjecture, and to extract various consequences from it. 

Eq. (1.6) incorporates our previous conjecture (1.5) when n = 1, as we discussed above. Using 
independent computations of a n [4, 18] and h n [16, 17], we have verified that (1.6) holds exactly for 
a free scalar and for a free massless Dirac fermion at integer values of n (up to n = 500) as well as 
in the n —> oo limit. Further, since we have a simple expression for h n that holds for any real n, we 
can predict the behavior of a n for non-integer n. This is plotted in Fig. 3 for the free complex scalar 
and Dirac fermion. Fig. 3 also includes cr n for a holographic CFT dual to Einstein gravity. This is 
a prediction of our conjecture as h n is known for these holographic theories, but direct access to cr n 
is currently limited to n — 1. Indeed, the celebrated Ryu-Takayanagi prescription [6] only provides 
a holographic description of entanglement entropy. 

A general formula for the scaling dimensions h n in dimensions d > 2 was obtained in [17] — see 

2 Eq. (2.5) provides the general equation relating d n h n \ n —i and C T in d dimensions. The above result follows from 
simply substituting d = 2 and also C T = c/(27r 2 ). 
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(2.4) below. As we explain in section 3, h n is expressed in terms of the thermal energy density £(T) 
of the CFT on the hyperbolic space Hk^ 1 (the simplest geometry with constant negative curvature). 
Using eq. (1.6), we find a new formula for the Renyi corner coefficient: 


a n = ^-R 3 (£(T 0 )-£(T 0 /n)) , 

71—1 


( 1 . 11 ) 


where R is the curvature radius of the hyperbolic plane HI 2 , which also determines the temperature 
To = 1/(27 tR). Note that neither cr n nor h n depend on R , and that this expression is always 
positive since the energy density £ (T) increases with temperature. Eq. (1.11) allows for explicit 
calculations of the Renyi corner coefficient a n for any n > 0, including non-integer values of n. ft 
also sheds light on the physical content of a n . For example, following the analysis of [17], it yields 
our earlier conjecture (1.5) where <J\ is proportional to C T , the central charge appearing in the 
vacuum two-point function of the stress tensor: 

{T^{x) T vk ( 0)) = %X^ K (x), (1.12) 

oc 

where X Alt/ ^ K (x) is a dimensionless tensor, whose structure is fixed by conformal symmetry [20]. As 
noted above, the desired result follows from expanding h n around n — 1 in eq. (1.6). Further, in 
the limit n —y 0, we find that the corner coefficient yields another important physical constant, the 
coefficient c s which controls the thermal entropy density of the CFT (in flat space), i.e., s flat = 
c s T d ~ l . In particular, we show that eq. (1.11) yields 


lim a n 

71 — 1-0 


c s 1 
127T 3 7l 2 


(1.13) 


In section 6 we extend this result by establishing that the corner function a n (8) diverges as 1/n 2 
for any opening angle. In the opposite limit, i.e., as n —> oo, eq. (1.11) becomes 


■Too = R 3 (£(T 0 ) - fill)) > 0 . 


(1.14) 


This quantity does not seem to lend itself to an interpretation in terms of flat space observables, 
in contrast to the n —> 0,1 limits. In section 6, we show that this asymptotic behavior extends to 
finite angles, i.e., a oc (9) is finite for all 6. In Fig. 2(b), we show the generic dependence of the a n 
coefficient on the Renyi index n as deduced from eq. (1.6). 

In Fig. 3, we have normalized the a n by the corresponding central charges C T , which makes them 
cross at n = 1 by eq. (1.5). The three curves all exhibit a similar dependence on n. In particular, 
they diverge as 1/n 2 for n —> 0, and asymptote to constant values for n —» oo. A salient feature of 
Fig. 3 is that the curves only coincide at n = 1, meaning that the ratio cr n /C T is only universal for 
the entanglement entropy, as was independently observed in [18]. 

Interestingly, the free fermion and holographic curves in Fig. 3 are hardly distinguishable in the 
whole range of values, with the agreement becoming even better for values 0 < n < 1. In particular, 
this shows that the ratios c s /C T are very close to each other in both theories. The exact values are 

C f C hDl 47T 5 

-4- = 12ttC(3)~ 45.3165, ~ 45.3362. (1.15) 

Om O7-1 Z I 
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Figure 3: Log-log plot of cr n /Cr versus the Renyi index n for the free complex scalar, free Dirac 
fermion, and holographic theories. Strikingly, the latter two are almost equal. We note that all 
three lines must cross at n = 1, where one recovers the universal ratio, oq/CV = vr 2 /24 ~ 0.411. 

While the irrational numbers involved in both cases are different, the final results differ by only 
approximately 0.04%. 

Our results for a n also reveal a surprising duality relating the corner coefficients of the complex 
scalar crjf and fermion od: 

= ( 1 - 16 ) 

which is valid for any real positive n. As we show in proving this result, it requires that the cor¬ 
responding thermal partition functions (on S 1 xB 2 ) are related in a relatively simple way. This 
surprising relation (1.16) also implies further connections between the two free theories, e.g., the 
above expression can be rewritten in terms of the conformal dimensions of the corresponding twist 
operators, yielding nh\ f ,f = — 

The remainder of the paper is organized as follows: In section 2, we review some results con¬ 
cerning twist operators and their conformal dimensions h n . In section 3, we compute h n for a 
free conformally coupled complex scalar and a massless Dirac fermion in three dimensions. We 
compare these results with those obtained previously for the Renyi corner coefficient <j n , finding 
perfect agreement with our conjecture (1.6). In this section, we also compute the coefficients K n , 
arising in the sharp corner limit (6 —> 0), for the free fields. These coefficients play a central role 
in section 6 — see below. Given our results for a n , in section 4, we discuss a surprising duality for 
the corner coefficients of the scalar and the fermion theories. In section 5, we use our conjecture 
to produce a new expression for a n in holographic CFTs dual to Einstein gravity. In section 6, 
we study the n dependence for corners of arbitrary opening angles. Then, for a given index n, we 
propose a method for estimating the curve a n {6 ) using the coefficients K n and a n alone. We apply 
this approach to construct interpolating functions for the scalar and the fermion for n = 1,2,3 and 
in the limit n —> oo, and we show that they accurately fit the lattice values in the cases in which 
these are available. A number of technical details are provided in the appendices: In appendix A, we 
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evaluate certain integral expressions for the scaling dimensions in the free field theories to express 
h n in terms of finite closed-form sums, valid for odd n. Using the integral expressions for h n , we 
also provide explicit values of a n for certain rational values of the Renyi index, which provides an 
interesting demonstration of the duality between the coefficients for the scalar and fermion theories. 
In appendix B, we describe the details of calculations for the K n presented in section 3. Appendix C 
explains the origin of a relation observed in [18] between the limit of the corner coefficient at large 
n and the thermal partition function on S l x HI 2 in the same limit. In appendix D, we provide a 
general proof of the duality between the corner coefficients of the free scalar and fermion theories. 
In appendix E, we connect the Renyi corner coefficients to the Renyi entropies S n for a circle, using 
the relation [17] between these entropies and the twist scaling dimensions h n . 

2 Twist operators 

As noted in the introduction, twist operators were originally defined in discussing Renyi entropies 
in two-dimensional CFTs [3]. Twist operators are well understood in this context since they are 
local primary operators. They can be formally defined for quantum field theories (QFTs) in any 
number of dimensions, with the replica method, e.g., see [16, 17, 21]. However, in higher dimensions, 
they become nonlocal surface operators and their properties are less well understood. The replica 
method begins by evaluating the reduced density matrix p v as a Euclidean path integral where 
independent boundary conditions are fixed on the region V as it is approached from above and 
below in Euclidean time, i.e., with t E —> (U. This path integral for the partition function Z n is then 
extended to the Euclidean path integral on a n-sheeted geometry [3] , where the consecutive sheets 
are sewn together on cuts running over V, to represent 

Tr W] = f|. (2-1) 

The denominator is introduced here to ensure the correct normalization, i.e., Tr[/?y] = 1. In 
defining the twist operator r„[U], the above construction is replaced by a path integral over n 
copies of the underlying QFT on a single copy of the geometry. The twist operator is then defined 
as the codimension-two surface operator extending over the entangling surface, i.e., the boundary 
of the region V, whose expectation value yields 

(r„)„ = TV[p?] . (2.2) 

where the subscript n on the expectation value on the left-hand side indicates that it is taken in 
the n-fold replicated QFT. Further, here and in the following, we omit the V dependence of r n to 
alleviate the notation. Hence eq. (2.2) implies that r n opens a branch cut over the region V which 
connects consecutive copies of the QFT in the ?7-fold replicated theory. Closely related to the twist 
operator, the so-called swap operator was introduced to compute Renyi entropies on the lattice 
[ 22 ]. 

In the case of a CFT, the conformal scaling dimension h n of the twist operator is defined by 
the coefficient of the leading power-law divergence in the correlator (T MJ/ r n ) n as the location of T fW 
approaches that of r n [16, 17]. In the case of a twist operator on an infinite (hyper)plane, as is 
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Figure 4: Planar twist operators r n in d = 2,3,4 used in the definition of the twist scaling 
dimension h n , via eq. (2.3). They have support on the boundary (in red) of the entangling region 
V, which is at time fg = 0. The stress tensor insertions lie a distance y from the twist operators. 
In d > 2, the use of planar boundaries here is a matter of convenience, as h n can be extracted from 
any entangling surface dV by taking y —>■ 0. 


illustrated in Fig. 4, this correlator reads 


(Tab 7~n)n 



_ d a b jrp \ p. 

X J J W ai Tn/n " , 

2tt y a 

h n (d — 1 )§ij — dhihj 
2n y d 


(2.3) 


where the indices i,j and a,b denote the two transverse directions and the d -2 parallel directions 
to the twist operator. 3 Further, y is the perpendicular distance from the stress tensor insertion to 
the twist operator and hi is the unit vector orthogonally directed from r n to the stress tensor. Note 
that here denotes the stress tensor for the entire n-fold replicated CFT. 

While the above expressions are only valid for a twist operator on a (hyper)plane, we stress that 
in general the leading singularity takes this form whenever y <C £, where i is any scale entering 
in the description of the geometry of the entangling surface. Hence the scaling dimension h n is a 
fixed coefficient which is characteristic of all twist operators r n (in a given CFT), independent of 
the details of the geometry of the corresponding entangling surface. Finally, let us add that hi = 0 
since the twist operator r n becomes trivial for n — 1. 


2.1 Insights from hyperbolic space 

In refs. [16, 23], the entanglement and Renyi entropies for a spherical entangling region V = S d ^ 1 
in general CFTs were studied using a conformal mapping from the conformal vacuum in flat space 
to a thermal ensemble on the hyperbolic cylinder. The Euclidean version of this transformation 
maps a n-fold cover of to the thermal spacetime S 1 x HI^ 1 . The flat space geometry has branch 
cuts beginning at a sphere of radius R , and in the hyperbolic geometry, the radius of curvature 
of Hk^ 1 is R and the period of S 1 is 2irRn. Using this construction, the conformal dimension h n 

3 Let us add that implicitly the above expressions are normalized by dividing by (r n ) n but we left this normalization 
implicit to avoid the clutter that would otherwise be created. 









can be expressed in terms of the thermal energy density £ (T) of the (non-replicated) CFT on the 
hyperbolic space HI f/ ~ 1 [16, 17]: 

/IT 77 

K = R d (£(To) - £(To/n )) , (2.4) 


where T 0 = l/(27ri£). Note that the conformal dimension is independent of R since this radius 
is the only scale in the above expression, i.e., both of the energy densities above are proportional 
to 1/R d . Of course, the R independence is required since h n is a dimensionless coefficient and as 
described above, it is independent of the geometry of the entangling surface. 

With the above definition (2.4), it is clear that setting n = 1 yields hi = 0. However, as noted 
in the introduction, studying the limit n —> 1 of this expression also allows one to prove [17] 

dnK Ui = + C T , (2.5) 


where C T is the central charged appearing in the two-point correlator (1.12) of the stress tensor. 

In the n —» 0 limit, the temperature diverges for the second energy density in eq. (2.4) and this 
contribution dominates the result for the conformal dimension. Further, at large temperatures, the 
curvature of the hyperbolic space becomes unimportant and to leading order this energy density 
will match that in flat space [24], In particular, for a d -dimensional CFT, the thermal energy 
density in flat space directly relates to the entropy density, s flat = c s T d_1 , which defines a useful 
coefficient c s that characterizes the number of degrees of freedom in the CFT. The corresponding 
thermal energy density in flat space is then given by £ flat (T) = c R^c s T d and hence, with £{T) ~ 
£ flat (T) [1 + 0(1/(RT) 2 )} at high temperatures, eq. (2.4) yields 


lim h n 

n—> 0 



( 2 . 6 ) 


The negative sign arises here because n < 1. Thus, in this limit, the conformal dimension of the 
twist operator diverges with (d- 1) power of the Renyi index, and is proportional to the flat-space 
thermal entropy coefficient. In particular, we can see this behaviour arises in two dimensions from 
eq. (1.9). The small-rt behaviour of the twist dimension is then 0 = — c/(12n) and comparing 
with eq. (2.6), we can extract c s = tic/3 , which precisely agrees with the well-known result for this 
coefficient in d — 2 [25]. 

In the opposite limit, n —> oo, the temperature vanishes in the second term in eq. (2.4), leaving 

/IT T) 

lim h n = —— R d (£(T 0 ) - £(0)) , (2.7) 

n—>-oc a — 1 

which shows that h n grows linearly with n when the Renyi index is large, independent of the 
dimension d. Let us comment here that the physical interpretation of the factor £ (T 0 ) — £(0) in 
eq. (2.7) depends on the details of the renormalization scheme used in defining the energy density 
£(T). For example, using a heat kernel regularization with free fields in the next section, we will 
hnd £(0) = 0, and hence eq. (2.7) yields h^ oc £(T 0 ). However, with the holographic calculations 
described in section 5, one finds £(T 0 ) = 0 and hence in eq. (2.7), one is left with hoc oc — £(0). 
In this case, £(0) < 0 is interpreted as the Casimir energy density for the boundary CFT on the 
hyperbolic space. In any event, these ambiguities for the values of the individual energy densities 
in eqs. (2.4) and (2.7) are eliminated in the difference £(T 0 ) — £(T 0 /n ) and so they do not affect 
the value of the scaling dimension of the twist operator. 
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3 Free fields 


In this section, we present explicit the free field calculations which verify our new conjecture (1.6) 
relating a n to h n . Hence in the following, we restrict the discussion to three spacetime dimensions, 
in which case, r n is a line operator as shown in Fig. 4. We begin by computing the scaling dimensions 
of the twist operators h n corresponding to a free complex 4 scalar and a free Dirac fermion in three 
dimensions. This allows us to simply compute the smooth corner coefficient a n for any real n > 0. 
We use this to confirm our conjecture in the case of free CFTs, and examine some physical limits 
of a n . In this section, we also compute the coefficients n n appearing in the opposite limit, i.e., that 
corresponding to a very sharp corner with 6 —)■ 0. We will use the latter in section 6 to estimate 
the value of the corner functions a n {6 ) for arbitrary values of the opening angle. 

3.1 Twist operators 

Recall that the scaling dimension of the twist operator r n in a general CFT can be evaluated with 
the expression in eq. (2.4), which in three dimensions reduces to 

h n = 7 thR 3 (. £(T 0 ) - £(T 0 /n)) , (3.1) 

where £(T) is the thermal energy density of the CFT on the hyperbolic cylinder S 1 x HI 2 . Recall 
that in this curved space geometry, R is the radius of curvature of HI 2 , and the period of S 1 is 
set to the inverse of corresponding temperature, with T 0 = 1/(27 tR). Now, our first step will be 
to determine the partition function of the theory on the hyperbolic cylinder and then the desired 
thermal energy densities can be evaluated as 

£{T) = t^' oeZ(T) ’ (3 ' 2) 

where is the regulated area of the hyperbolic plane HI 2 . To be precise, we have [16, 26] 

- 1 + 0{S/R)^j . (3.3) 

Through the conformal transformation discussed in section 2.1, UV divergences in the entanglement 
entropy or the Renyi entropy in M 3 become IR divergences of the corresponding (total) thermal 
entropy on the hyperbolic plane HI 2 . Therefore, in order to identify quantities in a meaningful way, 
the UV cut-off in the flat space calculations must be mapped to the IR cut-off in the calculations 
on the hyperbolic cylinder. Hence, we see the short distance cut-off 5 (in flat space) appearing in 
the regulated area above. 

The partition functions for free fields on the (Euclidean) manifold S 1 x HI 2 are readily calculated 
using heat kernel techniques and the results for the conformally coupled complex scalar and a 

4 Note that our discussion in [12] involved a real scalar, for which the values of h n and a n are half those given 
here. We consider a complex scalar here to facilitate the comparison with the free fermion. 


V H 2 = 27 rR 2 (j- 
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n 

n —> 0 

l+e 

2 

3 

4 

5 

6 

n —> oo 

K 

C(3) 1 

47 r 3 n 2 

128 

1 

24?r 

i 

27^/3 

37T+8 

1927T 

\J 25—2\/5 

125 

81+34\/37r 

1944tt 

Win 

y 

3t(3) 1 

7T f 

1 

5 

l+6\/2 

y/ 425+58\/5 

261+20\/3 

W 

167T 3 n 2 

128 

64 

108 \/3 

256 

500 

5184 


Table 1: Scaling dimensions h n of the twist operator at various values of n for a massless complex 
scalar and a massless Dirac fermion. In the third column |e| <C 1. For non-integer n, see Table 5. 


massless Dirac fermion are given respectively by [17] 5 


log Z CS (T) 
log Z\T) 


Vep T f°° du u 2 + sinh 2 u (u coth u — 2) 

2 R J 0 u 3 sinh 2 u smh(u/(2nRT)) 

VfipT r°° du 2 sinh u — u(u coth u + 1) 

2 R J 0 u 3 sinhw tanh(u/(27r.RT)) 


(3.4) 


where again R is the radius of curvature of the hyperbolic plane and the temperature T corresponds 
to the inverse period of the circle S 1 . It is now straightforward to evaluate the scaling dimension 
(2.4) of the corresponding twist operators using eq. (3.2) and the final expressions can be written 
as 


ic = 


du 


Ann 2 J 0 sinh u |_sinh 3 (w/n) 


cosh(w/n) 


— n 


h} = _ 

n 8vm 2 


du 


tanh u 


n 


2 + sinh 2 u 
sinh 3 u 


cosh u 
sinh 3 
2 + sinh 2 (w/n) 
sinh 3 (u/n) 


(3.5) 


In appendix A, we show that these integrals can be evaluated for odd values of n, and expressed in 
terms of sums involving trigonometric functions. In any event, the above expressions can be easily 
evaluated for different values of n, and we will use this to test the conjecture (1.6) relating h n to 
the corner coefficient a n in the next subsection. Some explicit results are collected in Table 1. 

Following our general discussion above, let us comment further on some of the results in the 
table. First, we consider the limit n —> 0 for which eq. (2.6) becomes 


lim h n 

n—> 0 


Cs J_ 
12n 2 n 2 


(3.6) 


By evaluating the integrals for small n, we find that our free field results have the desired form, 
with both of the scaling dimensions diverging as 1/n 2 — see Table 1. Further, examining the overall 
coefficient of these divergences, eq. (3.6) indicates that 


c 


CS 

s 


3<(3) 


and 


4 = 


9C(3) 

An 


(3.7) 


for the thermal entropy coefficient of the complex scalar and the Dirac fermion, respectively. One 
can easily verify that these values precisely match the expected thermal entropy coefficients for both 

5 We emphasize that both of these expressions are UV finite ( i.e the potential singularities as x —> 0 are removed). 
We also note that both yield £(0) = 0 and a nonzero value for £(Tq) — see comments below eq. (2.7). 
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fields, e.g., see [27]. Another interesting limit to consider is n —> 1. While the scaling dimension 
vanishes at precisely n = 1, evaluating eq. (3.5) in the vicinity of this point, we find 


h cs \ 

u n ln=l+e 
^ rl I n=l+e 


7r 



128 


17vr 2 

1920 6 
137T 2 

1920 6 


+ 0(e 3 ), 

+ 0(e 3 ). 


(3.8) 


Now recall that eq. (2.5) determined the leading coefficient in terms of the central charge C T . In 
particular, we have d n h n \ n= i = t^C t for d = 3. Now given C£ s ’ f = 3/(167 t 2 ) [20], this latter result 
yields d n h n | n=1 = 7r/128 for both fields, in agreement with the above expansions. Finally, the scaling 
dimension exhibits the expected linear growth with the Renyi index as n —> oo. In particular, we 
find: 


lim 1C 

n—> oo 


303 ) 

167T 3 


and 


lim h„ 

n—>oo 


03 ) 

47T 3 


n. 


(3.9) 


3.2 Corner Renyi entropies 

The regulator independent contributions to the entanglement and Renyi entropies produced by a 
sharp corner in the entangling surface for d = 3 free conformally coupled scalars and free Dirac 
fermions were studied in a series of papers by Casini and Huerta [4], However, the expressions 
for, e.g., the Renyi coefficients a n were left in a very complicated and implicit form. Motivated 
by our original conjecture (1.5) relating ny and C T , and the possibility that a similar universal 
relation could also exist for a n , the results for crjf ,f were recently reduced to simple and closed-form 
expressions at integer n [18]. These new expressions were used to easily evaluate the first coefficients 
with n > 1, as well as determine the n —» oo behaviour. In this section, we use these results to 
confirm our new conjecture, verifying that eq. (1.6) holds for both the scalar and the fermion. 

In the second part of this section, we consider the coefficients k™’, which control the corner 
contribution of the Renyi entropy in the limit 6 —> 0. In particular, we evaluate the corresponding 
expressions for n = 2,3,4, and at large n — see also appendix B. In section 6, we will use these 
coefficients, along with the crjf ,f , to construct simple interpolating functions to approximate the 
corner function aj) s,f (0) for all angles. 


3.2.1 Smooth surface limit 


The expressions for cr[f ,f appearing in [4] were recently shown to reduce to the following simple 
sums (for integer n) [18]: 


a 


CS 

n 



k =1 


k)(n — 2k) tan 
127rn 3 (n — 1) 


crl = 


(n—1)/2 

E 

k=—{n— 1)/2 


k(n 2 — 4k 2 ) tan (^) 
247m 3 (n — 1) 


(3.10) 


Now combining eq. (3.1) for the conformal dimension h n with our new conjecture (1.6), we find 
a new general formula for the corner coefficient in terms of the thermal energy density on the 
hyperbolic cylinder: 

77 

tfn =-- R 3 (S(T 0 ) - S(T 0 /n)) . (3.11) 

n — 1 
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n 

n —>■ 0 

l 

2 

3 

4 

5 

6 

n —y 00 

< 

C(3) 1 

i 

i 

i 

37T+8 

0 25—2\/5 

81+34\/37r 

3C(3) 

47t 4 n 2 

128 

247r 2 

54V3n 

576tt 2 

500tt 

9720tt 2 

167T 4 


3C(3) 1 

1 

i 

5 

l+6\/2 

1/425+58V5 

261+20\/3 

C(3) 

167T 4 n 2 

128 

647T 

216^3 ir 

768?r 

2000vr 

25920?r 

47T 4 


Table 2: Corner coefficient a n for various Renyi indices for the complex scalar and Dirac CFTs in 
d = 3. 


In particular, focusing on the free fields, we can use the expressions (3.5) for h n to produce the 
following prediction for the corner coefficients in these theories 

cs = i r du 

n 47 r 2 n 2 (n — 1) J 0 sinh u 


1 

r°° du 

0 2 + sinh 2 u 

71° 

2 + sinh 2 (u/n) 

87 r 2 n 2 (n — 1 ) 

) 0 tanh u 

1 L * Q 

sinh u 

sinh 3 (u/n) 


We have explicitly evaluated these integrals for a few values of n in Table 1, and also numerically 
plotted the coefficients for continuous n in Fig. 3. In appendix A, we also evaluate the integrals 
analytically for odd n > 1 to express these corner coefficients as certain closed sums — see eqs. (A.4) 
and (A.8). Our new closed-form sums agree with eq. (3.10) for all odd integer values that we have 
checked, 1 < n < 500, but take a very different form from the latter. It would be desirable to 
explicitly prove that both are equal for all odd integers. In any event, we have also verified that the 
above integrals (3.12) exactly reproduce the values obtained from eq. (3.10) for n — 1,2, ■ ,500, 

as well as in the limit n —> oo. In that limit, we have also verified that the subleading terms agree, 
as we discuss below. We find that this agreement provides strong evidence for our conjecture (1.6) 
relating the corner coefficient a n and the scaling dimension h n of the corresponding twist operators. 


cosh(u/n) 

sinh 3 (w/n) 


— n 


cosh u 
sinh 3 u 


(3.12) 


Of course, we can also use this agreement to combine the relation h n = ■n{n — l)cr n with eq. (3.10) 
to express the scaling dimensions in terms of sums: 


h 


CS 

n 



k= 1 


k){n — 2k) tan (^) 


hi 


(n-l)/2 

E 


k(n 2 


k=—(n— 1)/2 


4fc 2 )tan(^) 

24n 3 


(3.13) 


which is valid for all integer values of n > 1. 

One advantage of the integral expressions (3.12) for over the sums in eq. (3.10) is that 
they allow us to evaluate the corner coefficients for non-integer values of n. For example, we can 
examine how these coefficients vary in the vicinity of n — 1 , i.e., 


es I 

°n | n= l+ e 

f I 

°n ln=l+<E 


1 

128 

1 

128 


17 


1920 

13 


1920 


6 + 0(e 2 ), 

e + 0(e 2 ). 


(3.14) 


Hence we see that the agreement between the two coefficients at n = 1 does not extend beyond this 
precise point, as expected from eq. (3.8). We can also consider the limit n —> 0, which yields 


lim off ~ 

n-5-0 


«3)1 

47 t 4 n 2 ’ 


and 


lim al ~ 

n->0 


303 ) 1 

167 T 4 n 2 


(3.15) 
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Of course, these results are simply related to the corresponding limit of the scaling dimensions in 
eq. (3.6), and hence the corner coefficients are proportional to c s and diverge as 1/n 2 in this limit. 

Next, we note that while the h n diverge linearly with n as n —> oo in eq. (3.9), the a n asymptote 
to the constant values <r^ ,f given in Table 2. The detailed large-n asymptotic behavior of both 
corner coefficients can be determined from eq. (3.12) and reads: 6 


where 


lim <r„ 

rwoo 


<Go 





bf = 


u 2 


b? = -1 + 


7T 


45C(3) 


~ 0.8008, 


7 


12C(3) 


~ 0.3158, 


°4 




(3.16) 


(3.17) 


We further find that the series in brackets only contains even powers of 1/n, with terms that 
alternate in sign. These analytical results for the asymptotics agree with the numerical analysis 
presented in [18] of the sums in eq. (3.10). This analysis yielded precise estimates for b ™’, which 
coincide with the values given above. Based on their numerical fits, the authors of [18] conjectured 
that the terms in the large-n expansions appear in pairs with identical coefficients, which leads to 
the overall factor of 1 + 1/n in eq. (3.16). Here the large-n expansion of our integral expressions 
allow us to confirm the appearance of this factor analytically. As an aside, we observe that the same 
factor 1 + 1/n fully determines the n dependence of the universal coefficient of the d — 2 Renyi 
entropy in eq. (1.10). Overall, the large-n asymptotic analysis discussed above provides further 
support for our main conjecture (1.6) relating h n and a n . 

Ref. [18] further observed that for both of the free fields 


lim cr„ 

n—>oo 


1 

47T 2 


T 

J o< 


(3.18) 


where the coefficient T ac is related to the logarithm of the partition function on S 1 xH 2 , as calculated 
in [28]. Given the connection (1.6) between the corner coefficient and the scaling dimension, the 
origin of the above relation can be traced to the facts that the scaling dimension can be expressed 
in terms of the energy density on the hyperbolic geometry HI 2 , as in eq. (3.1), and that the energy 
density is determined by a derivative of logZ(T) with respect to the temperature, as in eq. (3.2). 
We defer the precise details to appendix C because producing the exact result in eq. (3.18) requires 
taking into account differences between our present conventions and those in [28]. Our analysis in 
the appendix establishes that eq. (3.18) is a general relation that holds for any three-dimensional 
CFT. 

6 We note that for the scalar, the large-n expansion can be performed directly on the corresponding integral in 
eq. (3.12). However, the analogous calculation for the fermion does not produce a consistent expansion. Rather in 
this case, we took advantage of an alternative representation of cr* which makes use of the expression for log Z f (T) 
given in [28] — see eq. (7.5) and appendix D. 
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3.2.2 Sharp corner limit 


The main goal of our paper is to establish the relation (1.6) between the scaling dimensions of 
twist operators and the corner coefficients cr n , arising for almost smooth corners. However, it is 
also interesting to study the opposite limit of the corner function a n (9), corresponding to a very 
small opening angle, i.e., 6 —y 0. In this limit, the corner contributions are characterized by the 
coefficients n n defined by 


a n (9 —t 0) — n n /6 . 


(3.19) 


We will evaluate the first few n ™' f here and use the results to produce simple functions to approxi¬ 
mate ajf ,f (0) for all angles, in section 6. Further, we will explicitly show that these simple functions 
provide an excellent approximation to the exact a n {6 ) for n — 1, 2, 3. 

Interestingly, in the case of free fields, the n n coefficients are related to the so-called “entropic 
c-functions” for a family of corresponding massive free fields in d = 2 [4], In particular, they can 
be evaluated as 


K cs,f _ 

7r 




where the ■ (t) are the universal functions 


C = 


(3.20) 


(3.21) 


where S n (rri£) are the Renyi entropies of d = 2 free fields with mass m for a single interval of length 
l. We leave the details of the calculations of the c^ s ’ f (t) to appendix B. The final expressions for 
the coefficients can be written as 


= 1 f°° dt t 2 u'k/n, cs - ^/n,cs(( 1 ~ 2fc / U f + ^ i 1 + U l/n, J) 

n i n - !) ^ Jo 2 ( 1 +«fc/n,cs) 

2 ( V /2 r,t ^ U hnM k / n f + ^ - U l/n, f)) 

~ <n-\) ^ J 0 2(1 — u 2 k/n i ) 


(3.22) 

(3.23) 


which applies for integer n > 1. The functions Uk/ n ,cs,f(t ) are solutions to the differential equations 
(B.4) and (B.6) given in the appendix, subject to the boundary conditions (B.5) and (B.7) in each 
case. Evaluating the expressions in eqs. (3.22) and (3.23) for n = 2, 3,4, oo, we obtain the numerical 
results shown in Table 3. 

As also described in the appendix, K cs,i corresponding to the special case n = 1 are evaluated 
using a similar approach using two-dimensional entropic c-functions. The values of these coefficients 
were evaluated by [4] and also appear in Table 3 under n — 1. 


3.3 n and 9 factorization? 

One can ask whether the Renyi index n and 9 dependences of the corner function a n (6) factorize, 
namely, 

a n (0) = /Hoi(6 ») , (3.24) 
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n 

l 

2 

3 

4 

n —> oo 

K CS 

0.0794 

0.0455996(1) 

0.037339(1) 

0.033798(1) 

0.0262(2) 

4 

0.0722 

0.0472338(1) 

0.040662(1) 

0.0376674(1) 

0.030(4) 


Table 3: Sharp-corner coefficient for various Renyi indices for the massless complex scalar and 
massless Dirac fermion in d = 3. 

where f(n = 1) = 1 and a\{9) is the corner function appearing at n = 1 (i.e., in the entanglement 
entropy). Such factorization occurs for the corner contribution in the extensive mutual information 
model [21, 29, 30], which can be written as 

a® xt (6 ) ) =- af xt (8) , where af xt (8) = C T [1 + (it — 6) cot 8} . (3.25) 

o n hn\n=i n — 1 8 

This result follows from a straightforward generalization of our calculations for n = 1 in [12] . We 
note that the precise form of h n is undetermined for this model, apart from the general constraints 
valid for any CFT. Factorization also occurs in a more concrete setting in a class of quantum critical 
Lifshitz theories with dynamical exponent z = 2, which have a corner function characterizing their 
corner entanglement just like CFTs. Moore and Fradkin [31] indeed found that f(n ) = 1, i.e., the 
corner function a n {9) is independent of n for all angles 7, 8 . 

Using our above results, we find that the factorization (3.24) does not hold for the free complex 
scalar or for the Dirac fermion. Indeed, if eq. (3.24) were to hold, we would have the equality 
K n / k i = (7 n j(j\. It is sufficient to examine the n = 2 case to fold a contradiction: 

«£■/«? = 0.574 , <rf/<T? = 0.540 , 

/4M = 0.654, a l 2 /a[ = 0.637. (3.26) 

The deviations from factorization here are relatively small, i.e., about 6% for the complex scalar and 
3% for the fermion. However, these deviations are much greater than the numerical uncertainties 
in our evaluation of the corresponding K n coefficients — recall that the a n are exact. 

The failure of the factorization here is reminiscent of a similar non-factorization found in [16] for 
the Renyi entropies of spherical regions in higher dimensional CFTs. There the authors considered 
the question of whether or not the coefficient of the universal contribution in these Renyi entropies 
took the form f d (n ) a d , where a* d is a constant determined by the specific data of the underlying 
CFT but fd{ri) is some function of the Renyi index which is the same for any CFT in d dimensions. 
This kind of factorization is found in the Renyi entropy of a single interval in two dimensions, as 
can be inferred from eqs. (1.9) and (1.10). However, by examining a variety of holographic models, 
ref. [16] found that this factorization fails in any number of spacetime dimensions beyond two. 

7 We thank S. Furukawa for pointing this out. 

8 It has been numerically shown that non-analycities can arise in the n-dependence of the corner function of certain 
Lifshitz QCPs when these are studied using lattice realizations [32, 33]. However, there is no evidence that such 
singular effects arise in the ideal (fixed point) held theories describing those QCPs [31]. We thank Joel Moore for 
bringing this point to our attention. 
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4 Bose-Fermi duality 


In this section, we study duality relations between the universal entanglement coefficients at Renyi 
indices n and 1/n, in both d — 2 and 3, the latter being our main focus. First, let us make some 
observations comparing the corner coefficients cr n of the complex scalar and the Dirac fermion in 
d = 3. We note that the ratio ujf/cx^ monotonically decreases as n runs from 0 to oo, i.e., 


a cs 4 

°n-> 0 _ ^ 
f O > 

a n-s>0 6 


or _ 

j i 


<7 


1 


OO 

T~ 


3 

4 ’ 


(4.1) 


Interestingly, the values corresponding to the complex scalar and the Dirac fermion are closely 
related in the opposite limits of small and large n: 


n 2 a c J I 


>o 700 ’ 


n 2 cr f 


n I n—>0 


= CD 


(4.2) 


This is not coincidental as these relations constitute special cases of a new duality between the 
corner coefficients of the two free CFTs: 


n 2 a™’ = a i /n> n > 0, n G M. (4.3) 

The duality interchanges the boson and fermion coefficients, together with n 1/n. In particular 
crjf is fully determined by a/, and vice versa. In appendix D, we prove this relation for any positive 
real Renyi index, assuming our conjecture (1.6). As an immediate consequence of this duality, we 
can use our large-n expansion eq. (3.16) to get the following expansion in the opposite limit: 


cs, f 
°n —¥0 ~ 


Ct 


f, CS 


n* 


(1 + n) 


-I | rf.cs 2 

1 + o 9 ’ n 


'}a cs n 4 


+ 0(n 6 ) 


(4.4) 


with the same coefficients &j) cs as in eq. (3.17). Of course, combined with eq. (3.6), these results 
imply that the coefficient of one theory is determined by the thermal entropy coefficient c s of 
the “dual theory”. 

This surprising duality between the corner coefficients in the two free theories has numerous 
other consequences. For example, we can rewrite it in terms of the conformal dimensions of the 
corresponding twist operators as 


n h 


CS, f 

n 



(4.5) 


We can go further and use eqs. (1.13) and (1.14) to establish relations between the thermal entropy 
coefficient of each field with the thermal energy density on the hyperbolic plane HI 2 at To of the 
other: 9 


of, cs 


£ c ' (To) = 


12vr 3 R 3 ’ 


(4.6) 


where again To = l/(27rR). In fact, combining eqs. (3.1) and (4.5), we find a more general relation 
between the energy densities of the two theories at generic temperatures on the hyperbolic cylinder, 


T cs ’ f (T 0 2 /T) - £ cs,f (T 0 ) = (£ f,cs (T) - T f ’ cs (T 0 )) . (4.7) 

1 o 

9 The form of this result relies on £ cs ’ l (T = 0) = 0 in our conventions. In appendix D, we have adopted the 
conventions of [28] in which f cs,f (T 0 ) = 0 instead. As explained below eq. (2.7) and in the appendix C, the 
differences arise from the choice of the renormalization scheme for the partition function on S 1 x H 2 . 
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In appendix D, we establish a general relation for the logarithm of the partition functions of the two 
theories on S 1 x HI 2 and so there may be interesting relations, similar to eq. (4.7), for other thermal 
quantities. 111 These relations are relating the properties of one theory at temperature T to the other 
at Tq/T. Hence, in this regard, our duality is somewhat analogous to the Kramers-Wannier duality 
[34], which relates quantities in the two-dimensional Ising model at low and high temperatures. 

We point out that a duality closer in spirit to Kramers-Wannier can be obtained by considering 
the free CFT containing both the complex scalar and the Dirac fermion. This CFT has AT = 2 
supersymmetry, as it consists of a free chiral multiplet. Now, the corner function a n {6) is additive 
in that case as free fields factorize, implying a ^ ee ^= 2 = cbf + cd. The above Bose-Fermi duality (4.3) 
thus leads to 


n 2 a l ™ M=2 = > n > 0, n e M , (4.8) 

which is a self-duality connecting a lreeAfc2 at q = n and 1 jn. In particular, in this theory, 
acquires a concrete physical meaning as it is determined by the flat-space thermal entropy coefficient 
c s of the same theory. It would be interesting to investigate the fate of this self-duality in the 
interacting J\f = 2 supersymmetric CFT, i.e., the chiral Wess-Zumino model, which can be obtained 
by perturbing the theory of the free chiral multiplet. 

It is further natural to ask: Do similar entanglement dualities exist in dimensions other than 
d — 3? Below we point out that this indeed occurs in d — 2. 


4.1 Entanglement self-duality in d = 2 

For any CFT in d = 2, we have the following self-duality relation for the twist dimension: 

K' = ~Kjn • (4-9) 

as follows from eq. (1.9). This can then be rewritten in terms of the coefficient <r[ l 2) of the logarithmic 
contribution to the Renyi entropy of a single interval as 


( 2 ) - ^ 


na y = a 


1/n 


where we recall from eq. (1.10) that 


/ii 2) 


(2) = --- = ( 1 + Z. 

n n — 1 12 1 


(4.10) 


(4.11) 


The relation (4.10) is clearly reminiscent of the d = 3 Bose-Fermi duality presented above, and in 
particular of its manifestation as the self-duality in the free supersymmetric theory in eq. (4.8). As 
a special case of (4.10), we have a Bose-Fermi duality between a free Dirac fermion and a real scalar 
(compact or not): na ( ^ ),s = crW , as both have c = 1. 

Given the above results in d = 2,3, we are led to propose that relations of the form rA~ 1 o ( ^ ] ~ 
(j[y n or n d ~ 2 h 4 d) ~ — might occur as well in other higher-dimensional theories. The dualities 
may connect n to 1/??. in two different theories, as in eq. (4.3), or within the same theory, as in 
eqs. (4.8) and (4.10). 

10 There may also be interesting relations for the circle Renyi entropies using the results in appendix E — a topic 
we leave for future work. 
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n 

n —> 0 

1 

2 

3 

4 

5 

6 

n —> oo 

^.hol //^ihol 

a n /u T 

7r 2 

7r 2 

7r 2 (5\/l3—1) 

tt 2 (25V7-1) 

5tt 2 

7r 2 (73v/l9—1) 

7r 2 (53Ul09-l) 

7r 2 

81 n 2 

24 

648 

2916 

243 

16200 

29160 

36\/3 


Table 4: Corner coefficient for various Renyi indices in holographic CFTs dual to Einstein gravity. 


5 Holography 

Let us now consider strongly coupled holographic CFTs dual to Einstein gravity. We will use 
known results for the twist dimension h^° l to obtain the n-dependence of the corner coefficient a]) 01 . 
We begin with the AdS/CFT correspondence in its simplest setting, where it describes a given 
d-dimensional boundary CFT in terms of (d + l)-dimensional gravity in the bulk with the action 

'd(d- 1) 


I = 


1 

167tG 


d d+1 x^~g 


L 2 


+ n 


(5.1) 


where G is the (d + l)-dimensional Newton’s constant, L is the AdS^+i radius, and 1Z is the Ricci 
scalar. In order to obtain h^° l , we need to consider the thermal ensemble of the CFT on the 
hyperbolic geometry appearing in the construction of [23], which is then equivalent to a topological 
black hole [35] with a hyperbolic horizon. Using this connection together with eq. (2.4), the scaling 
dimension of twist operators in the corresponding holographic CFTs was evaluated as [16] 


hl° l = 


d -1 


8 G 


n 


y/d 2 n 2 — 2 dn 2 + 1 + 1 
dm 


d -2 


( Vd 2 n 2 — 2dn 2 + 1 + l)‘ 
d 2 n 2 


(5.2) 


In the following, we will focus our attention on the case of a three-dimensional boundary theory. 
Of course, one particular example of such a boundary CFT would be the supersymmetric gauge 
theory constructed in [36], in an appropriate large N and strong coupling limit. 

Setting d = 3 in eq. (5.2) and using the conjectured relation (1.6), we can predict the following 
form for the Renyi corner coefficient in these holographic CFTs: 

, 2 (3n 2 — 2)\/l + 3n 2 — 2 


u ho1 = 


108ttG 


n 2 (n 


i) 


(5.3) 


Of course, at present, there is no method available by which we may evaluate these coefficients 
directly in a holographic framework. This expression can be easily used to compute the corner 
coefficient for arbitrary values of n, and some results are presented in Table 4. In the table, we 
normalized the corner coefficient by dividing by the central charge C^ oX = 3L 2 /(tt 3 G) [37, 38]. We 
recover the expected result for the corner coefficient in the entanglement entropy (he., with n — 1). 
That is, a\ o[ /C'^ o[ = 7t 2 /24 [12, 13]. Let us note at this point that the numerical values for the 
<r n /C T in the holographic CFT are remarkably similar those for the free Dirac fermion, as shown in 
Fig. 3. The relative error is no more than 2.6% for n > 1 and 0.2% in the range 0 < n < 1. 

Turning now to the limits n —>• 0 and n —)• 1, we find 


h°i _ 
°n-s>0 ~ 


cr. 


hoi I 


I n= 1+e 


27t tG 
L 2 
8 ^G 


4 + - + <?( 

7i- n 

1 - le + 0(e 2 
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') 


(5.4) 
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Comparing the n —> 0 limit above with eq. (1.13), we find 


hoi 

L s 


4tt 2 L 2 
9 G 


(5.5) 


This result exactly matches the thermal entropy coefficient for the holographic CFT calculated from 
a planar AdS 4 black hole in the bulk, e.g., see [13]. 

In the large-n limit, we find: 


where 


hoi 


= cr 


hoi 

OO 




2a 


hoi 

OO 


3y/3n 3 ’ 


i>2 



—- - - ~ 0.01 and cr^ 1 
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We emphasize that we have shown analytically that the overall factor of (1 + 1/n) appears in the 
first term of eq. (5.6). Of course, this was the structure conjectured for the large-n expansions in 
eq. (3.16) for the free scalar and fermion CFTs. However in contrast to the free CFTs, the second 
term proportional to 1/n 3 spoils this simple n-dependence for the holographic CFTs. As explained 
in eq. (1.14), the coefficient a^ is proportional to S(T 0 ) — £( 0), a difference of thermal energy 
densities on the hyperbolic cylinder. We note that the holographic calculations yield £(T 0 ) = 0 and 
£(0) = — L 2 /(12v / 3vrGi? 3 ) here [39]. These results contrast to those for the corresponding free held 
energy densities in section 3, where £(0) vanishes while £(T 0 ) does not. 

Let us close this section by mentioning that our results here should be easily generalizable to 
other holographic theories. In particular, it would be interesting to study these corner contributions 
for holographic theories whose bulk action contains higher-curvature terms, e.g., using the results 
of [16]. 


6 Corners with arbitrary opening angles 

In sections 3 and 5, we saw that the corner coefficient a n has an interesting systematic dependence 
on the Renyi index. Below, we argue that this dependence is a general property of the entire 
corner function a n (9), at all opening angles 9. Further, in this section, we provide an efficient 
approach to approximating the angular dependence of the corner function a n (9), using knowledge 
of the coefficients tz n and a n , which control the asymptotic behaviour in the limits 9 0 and 7r, 

respectively. 

6.1 General n dependence 

We now go beyond the nearly smooth limit and establish some general properties of the Renyi index 
dependence of the corner function a n {9 ), valid for any opening angle 9. Without loss of generality, 
we consider a region V containing a single corner with opening angle 0 < 9 < tt but that is otherwise 
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smooth, as in Fig. 1(a). As it will be useful below, we recall that the Renyi entropy associated with 
V is then given by an expression of the form shown in eq. (1.2), i.e., 

Sn(V) = B n - a n (9) log (A(dV)/8) + • • • , (6.1) 

o 

where we have written the leading term in a geometric fashion, with A(dV) being the total length 
of the entangling surface, and 5 a UV cutoff. 

Let us begin by considering the n —> 0 limit of the corner function. Here we begin by adapting 
the arguments that were encountered in section 2.1 for the present discussion. First, we step 
back to consider a circular region with an entirely smooth boundary. The problem of evaluating 
the Renyi entropy for such a region in the vacuum of a general CFT can be mapped to one of 
examining a thermal ensemble on hyperbolic space HI 2 . However, it was observed by [24] that in 
the n —>■ 0 limit the result is dominated by the high temperature behaviour, which in fact matches 
that in flat space. Hence we have that the Renyi entropy for a circular entangling curve in three 
dimensions diverges as S n ~ 1 /n 2 . In particular then, this argument indicates that as n —» 0, the 
coefficient B n in eq. (6.1) diverges as 1/n 2 for a circular entangling surface. However, this coefficient 
is independent of the details of the geometry and so the same divergence appears in the leading area 
law contribution to the Renyi entropy for any entangling surface. Now note that the corner makes 
a negative contribution to S n with the second term in eq. (6.1) but the total Renyi entropy must 
always be positive. Hence we conclude that as n —* 0, a n (9) can diverge at most as fast as l/n p , 
with p < 2. Indeed, as 9 —> tt, we have a n (9 ) ~ a n (9 — 7r) 2 where a n diverges precisely as 1/n 2 , 
eq. (1.13). Next, we can use the fact that a n {9 ) decreases monotonically for 0 < 9 < tt, eq. (1.3), 
to establish that for all angles a n {9) diverges at least as fast as 1 jn v with p > 2. Combining these 
two inequalities we are lead to p = 2, i.e., the corner function diverges as 1/n 2 for any 9 as n —> 0. 

Next we consider the opposite limit, n —» oo. In this case, we know that a n asymptotes to a 
finite value a ^ > 0. Hence at this point, we again invoke the monotonically decreasing behaviour 
of the corner function to infer that a n ^. 00 (9) must either tend to a finite value or diverge, i.e., it will 
behave as n q with q > 0 in the large-n limit. However, if a n ^. 00 (9) was growing (with q > 0), the 
positivity of S n would require that the coefficient B n of the area law term and hence the entire Renyi 
entropy are also growing with n at least as fast as the corner function. However, it follows from 
considerations of information theory that S n must satisfy various inequalities [40] and in particular, 
one finds that d n S n < 0. That is, the Renyi entropy can not grow with increasing n (at any point). 
Therefore, we are led to the conclusion that the corner function remains finite as n —>■ oo (i.e., 
q = 0), again for any opening angle. Our calculations for the free CFTs confirm this result since 
the sharp corner coefficients K n of the complex scalar and the Dirac fermion approach finite values 
as n —» oo, see section 3.2.2. 

The general asymptotics of the corner function found above for the limits n —> 0 and oo, as well 
as the behaviour for large and small angles, are schematically illustrated in Fig. 2(a). 

6.2 Simple approximations for a n (0) 

In this subsection we address the following question: For a given Renyi index n, can one efficiently 
approximate the angle dependence of the corner function a n (9)l We will show that by knowing 
the asymptotic coefficients n n and cr n alone, one can construct a simple closed-form expression for 
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a n {6 ) that is not only very accurate for all 9 , but also exact as 9 —» 0, n. We will use two physically 
motivated functions that obey the correct asymptotics as our building blocks: 


a w (9) 


( Q -*) 2 

9(2,71 - 9) ’ 


cP\9) = 1 + ( 7 t — 9) cot 9. 


( 6 . 2 ) 

(6.3) 


The first one is the simplest algebraic function with a quadratic zero at 9 = n and with the pole 
at 9 = 0, as a n (9 ) is expected to have on general grounds — see eq. (1.4). The additional factor 
of (27T — 9) in the denominator ensures the function satisfies the constraint 0,^(9) = — 9), 

as required for the Renyi entropy of pure states. Interestingly, precisely the same function yields 
the corner function a n (9 ) of certain Lifshitz quantum critical points [31]. The second function a ^ 
corresponds to the functional form of the corner contribution in the extensive mutual information 
model [21, 29, 30], as described above at eq. (3.25). Both functions satisfy the expected constraints in 
eq. (1.3). However, we note that they also satisfy a stronger constraint obtained for the entanglement 
entropy coefficient di{9) using strong subadditivity and Lorentz invariance [4, 41], i.e., dgd^\9) > 
|<9e>a^(0)|/sind for 0 < 9 < tt. Further, as mentioned in the introduction, reflection positivity 
also imposes a infinite tower of nonlinear constraints on the corner functions a n {9) [10]. These 
constraints can be schematically written as det({c>^ +A+ 2 a n ( 6 l )}^,'l 1 0 ) > 0, with M > 1 and 9 € [0, 7 r], 
For example, with M — 1, this reduces to the linear constraint dgd n > 0, but the second inequality 
with M = 2 reads: 


dga n dgd n ~ (d^dn) 2 > 0 . (6.4) 

We have managed to explicitly verify that our basis functions d^\9) and d ( - 2 \9) in eq. (6.2) satisfy 
these non-linear constraints for M — 2,... ,6. 

We now construct an approximation d n (9) for the corner functions with a simple linear combi¬ 
nation of and a*- 2 -*: 


a n (9) ^ a n (9 ) = A.[ l 1) a (1) (6») + A^V 2 ^), 


(6.5) 


where the constants A $ are fixed by requiring the right-hand side of eq. (6.5) to have the same 
asymptotics as d n (9 ) in the limits 9 0 and tt. Taylor expanding a>( 1 \d ( ' 2 ' 1 in the these limits, we 

obtain: 


0 —t 0 : — A^ + 7 tA^ — K n , 

9 —t tt : —j^A^ + -A® = a n . 

7T Z 3 

The unique solution of eq. (6.6) for the two unknown coefficients is then: 

\(1) O ^ n 37T(J n \(2) ^ 2,H n TT <J n 

^n ~2 c J ^n ~Z ~n 

7T Z — 6 7T 7T Z — 6 


( 6 . 6 ) 


(6.7) 


The holographic CFT dual to Einstein gravity provides a simple test of this approach (6.5), at 
least for n — 1. For this theory, we have an integral expression to evaluate a) 101 ^) for any angle 
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Figure 5: Holographic corner function a\ ol (9) [42] normalized by C^ 61 , and the approximate form 
obtained by interpolating the asymptotics with eq. (6.5). The two lines can hardly be distinguished. 
The inset shows the ratio of the exact and extrapolated functions. 




(a) (b) 


Figure 6: Closed-form approximations for the corner function a n (6) at n = 1,2,3 and n = oo for 
(a) the free complex scalar and (b) the free Dirac fermion. The dots were obtained from numerical 
lattice calculations [4], The approximate functions obey the exact asymptotics as 6 —> 0,7r. 
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[42], and the necessary coefficients are n\° l = T(3/4) 4 L 2 /(2ttG) and cr]) 01 = L 2 /(8nG), e.g., see [13]. 
In Fig. 5, we show the curves for the exact result for a\ ol {9) / C^° l and the approximate function also 
normalized by C^° x = 3L 2 /(n 3 G). The agreement is seen to be excellent, with the relative error 
being less than 1% over the entire range of angles. 

We can also use the values of a n and n n computed in the section 3.2 to construct approximations 
of a n {9) for the free complex scalar and the free Dirac fermion. Now in principle, the a n {6) functions 
can also be determined exactly for these free fields, as described in [4]. Unfortunately, their results 
are written implicitly in terms of several functions which must be determined through a complicated 
system of coupled non-linear algebraic and differential equations. We explicitly used this approach to 
test our approximation for a!f (0) for the complex scalar, and we found that the difference between 
our approximation and the ‘exact’ curve was comparable to our numerical errors in evaluating 
a n {6) with the method of [4], The numerical errors were of the order of 0.20% while the maximum 
difference between the two curves was approximately 0.33%. Clearly, these results suggest that 
our approach produces a faithful approximation to a n [6) for the free fields. In Fig. 6, we show our 
approximate corner functions for n = 1, 2, 3 and oo, both for the free scalar and fermion. The figure 
also includes the results of lattice simulations for 6 = 7t/4, 7t/ 2 and 37 t/ 4, in the cases where these 
are known, 11 and again, we see that our interpolating curves agree very well with these lattice points. 
As an additional check, a c l(-K/2) was recently estimated to be 0.0086(2) by numerically studying the 
Goldstone phase of antiferromagnets on the square lattice [9]. Using a different system, a numerical 
linked cluster expansion [43] yielded 0.0094(2). These values are in excellent agreement with our 
estimate 0.0094, obtained with the simple interpolating method described in this section. 

Finally, we comment on the reasons behind the excellent accuracy of our interpolations. As 
noted in eq. (1.3), a n {6) is a monotonic convex function of 6. Thus we may expect that once the 
behavior in the vicinity of the endpoints 9 = 0 and 7r is fixed (by our knowledge of K n and <r n ), 
little freedom remains at intermediate angles. Our findings confirm this suspicion. Besides, we have 
checked that all the interpolating functions shown in Figures 5 and 6 satisfy the first six inequalities 
arising from reflection positivity mentioned above. 


7 Discussion 

In this paper, we have introduced a generalized version of the conjecture presented in [12, 13] for 
the coefficient appearing in the universal contribution in the entanglement entropy for an almost 
smooth corner, 

(7 - 1) 

In particular, we proposed that the analogous Renyi corner coefficients a n — see eq. (1.4) — are 
related to the scaling dimension h n of the corresponding twist operators r n by 

1 h n 

7T n — 1 

11 We thank Horacio Casini for providing these lattice results, as well as his assistance in directly evaluating a§ s (9) 
using the method of [4]. We obtained a™(0) for 0.6 < 9 < 3.0, but do not show it in Fig. 6 since given the uncertainties 
described in the main text, the result would be indistinguishable from the approximate interpolation in the plot. 
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This expression reduces to our previous conjecture in the limit n —» 1, and is therefore supported in 
this limit by the extensive results presented in [12, 13, 14, 15, 18], corresponding to free scalars and 
fermions, as well as a general proof for holographic theories. Here we have checked the generalized 
version of the conjecture for a free conformally coupled complex scalar and a free massless Dirac 
fermion for integer Renyi indices 1 < n < 500, as well as in the limit n —)■ oo. The calculations 
involved in the evaluation of a n and h n for the free fields are very different, both conceptually and 
computationally, which motivates us to conjecture that eq. (7.2) is a universal expression for general 
three-dimensional CFTs. It would also be interesting to explore how these expressions get modified 
in non-conformal theories, see e.g., [31, 44, 45]. 

Further, based on the properties of h n at small n, we have obtained that 


lim <r n 

71 — 1-0 


1 Cs 

127T 3 n 2 ’ 


(7.3) 


which yields a key observable, namely the thermal entropy density coefficient c s of the original CFT. 
As an interesting application of this result, we find that a n for the 0(1V) Wilson-Fisher IR fixed 
point in the N —)■ oo limit differs from its value at the UV Gaussian fixed point. Indeed, the ratio 
of the entropy coefficients at N —» oo is Cg R /c RV = 4/5 [46]. Hence a n already distinguishes the 
entanglement structure of these two CFTs even in the absence of finite N corrections, in contrast to 
the RG monotone F [47] . 1J This suggests that the full corner function a n {6) differs at other values 
of n and 6 between the A r —)-oo IR and UV fixed points. It will be interesting to verify this claim. 

Our original motivation for the new conjecture (1.6) was as follows: Essentially, the corner 
coefficient a n describes the universal response of the Renyi entropy to a small deformation of the 
(originally smooth) entangling surface, as shown in Fig. 1(b). This small displacement of the twist 
operator can be accomplished by making appropriate insertions of the stress tensor, following the 
approach of [48] — see also [49, 50]. Hence in the case of the entanglement entropy, the response 
is determined by correlators of the stress tensor and the entanglement Hamiltonian [48, 51]. Since 
the latter is given by an integral of the stress tensor over the region V [23] , the calculation actually 
involves correlators of the stress tensor with itself and hence it is natural that the response is 
controlled by the central charge C T as in eq. (1.5). In the case of the Renyi entropy (with n / 1), 
the corresponding calculation instead involves correlators of the stress tensor and the twist operator 
and so this reasoning naturally suggests that the response would be determined by the scaling 
dimension h n as in eq. (1.6). This reasoning also outlines a path potentially leading to a proof of 
our new conjecture, however, the detailed calculations are subtle and we must leave this topic for 
future work [52]. 

In eq. (1.10), we observed that the universal contribution in the Renyi entropy of an interval 
in a two-dimensional CFT also contains a factor of h n /(n — 1). This suggests that the form given 
in eq. (7.2) generalizes to higher dimension as well. In particular, using the results of [50, 53], we 
have shown that the higher dimensional analogs of o\ describing the universal contribution in the 
entanglement entropy for an almost smooth cone is proportional to the central charge C T for general 
holographic CFTs in any number of dimensions [15]. 13 Further, as shown in eq. (2.5), d n h n | re=1 is 
also related to C T in general dimensions. Hence, it is natural to suggest that eq. (7.2) will have a 
higher dimensional counterpart for the coefficients cx/ d) controlling the response of the Renyi entropy 

12 Of course, the scaling dimension h n of the twist operator serves the same purpose. 

13 Similar observations were recently made in [14, 44]. 
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to an almost smooth conical defect in the entangling surface of the form 


°n' 1 = 9(d) 


K d) 

n — 1 ’ 


(7.4) 


where g(d) is a geometric factor that depends on the spacetime dimension, and h^ ] is the scaling 
dimension of the corresponding twist operators in d dimensions. Proving this conjecture would give 
a new perspective on the geometric character of Renyi entropies in higher dimensions [49, 50, 54]. 

In section 4, we observed an intriguing duality between the corner coefficients of the free scalar 
and the free fermion. In particular, this relates the corner coefficients for inverse values of the 
Renyi index through the simple expression n 2 chf ,f = cr^, which we proved in appendix D. Further 
implications of this duality include that the energy density of one theory on the hyperbolic cylinder 
at temperature T is related to the corresponding energy density of the dual theory at T 0 2 /T, as in 
eq. (4.7). It would be interesting to explore the physical implications of this new duality in more 
depth. Similarly, it would be worthwhile to search for similar relations for other theories in other 
spacetime dimensions. 

Our proof of the Bose-Fcrmi duality makes use of the representation of log Z CS} f given in [28] - 
see appendices C and D. This allows us to give the following alternative representation of cr“ ,f 


cr 


CS 
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n 
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47r 2 (n — 1) 
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47r 2 (n — 1) 
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47T 2 


2 



d\ A 2 tanhA (coth(nA) — 1) , 
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7T 2 7T 2 



(7.5) 


It would be interesting prove the equivalence of these integral expressions with those in eq. (3.12), 
which were used in the main text. Alternatively, it also appears that these expressions may be 
somewhat simpler to relate to the sums in eq. (3.10), than those in eq. (3.12). 

In this section 6, we devised an approach to constructing a simple analytic function which gives 
a (very) good approximation to the entire corner function a n (8) for all angles. These functions make 
an informed interpolation between the asymptotic limits (1.4) at 9 —» 0 and 7r. Hence the only 
data necessary to construct these approximations is the almost smooth corner coefficient o n and the 
small angle coefficient K n in eq. (1.4). Of course, our conjecture (7.2) indicates that the coefficient 
o n can be accessed by a different calculation. In particular, the scaling dimension h n is determined 
by the leading singularity in correlation function of the stress tensor with a twist operator, as shown 
in eq. (2.3). Similarly, the small angle coefficient K n can be accessed by an alternative calculation. 
In particular, as described in [4, 13], the same coefficient controls the universal contribution in the 
Renyi entropy of an narrow strip. It would be interesting to investigate situations where these two 
alternate calculations would allow the corner contribution to be determined more simply than by 
evaluating the Renyi entropy for entangling surfaces with corners of varying angles. 


Note added — After the first version of this paper appeared, Dowker [55] presented an analytic 
proof of our conjecture (7.2) in the special case of a free scalar using eq. (3.10) and his own results 
for the scaling dimensions h r n s . 
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A Evaluation of h n for the free CFTs 

In this appendix, we evaluate the integral expressions in eq. (3.5) for the scaling dimension of the 
twist operator with odd Renyi index n > 1, in the free conformal field theories. Using the same 
integrals along with eq. (1.6), we also evaluate a n in closed-form for various rational non-integer 
values of n. These results provide an interesting explicit check of the Bose-Fermi duality discussed 
in section 4, and proven in appendix D. 


A.l Complex scalar 


We begin with the conformal dimension h n in the theory of a free complex scalar, which may be 
written as: 


h 


CS 

n 


l 

87m 2 



du 

sinli u 


cosh (-u/n) q coshu 
sinh' {u/n) sinh u 


(A.l) 


Here we have used the fact that the integrand is even to extend the range of integration over the 
entire real line. In passing, we also note that the integrand is finite at the origin, where it becomes 
(n 4 — l)/(15n) + 0{u 2 ). Further we note that the integrand decays exponentially fast as u —> ± 00 . 
Next we observe that if we shift the contour in the complex u plane from the real axis to A(n) = inn, 
then the integrand is invariant up to an overall phase of e m ' K = (—) n . Hence for odd n, we may 
re-express the conformal weight in terms of complex contour integral which runs along the contour 
C , illustrated in Fig. 7: 
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CS 
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1 

167m 2 



du 

sinh u 


cosh (u/n) , cosh u 

-77——— — n 6 -o— 

sinh ( u / n ) sinh u 


(A.2) 


Thus the result will be expressed in terms of the residues of the poles found to lie within the 
contour C. The integrand within the contour is only singular along the imaginary axis where sinh u 
vanishes, i.e., at u = ink with k — 1, 2, ■ • • , n — 1. The second term in the integrand actually has 
a quartic and quadratic pole at each of these points but neither of these contribute to the closed 
contour integral. Hence only the first term contributes with the result 
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CS 

n 


1 yv u-+i co s(nk/n) 
8n 2 sin 3 (7r k/n) 


(A.3) 
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Figure 7: The integration contour C in the complex u plane. Poles occur at u — ink, for all 
integers k that are not multiples of n. Here we illustrate the case n = 3, with two poles lying inside 
C. 


We note again that this result only applies for odd values of n. One can confirm that this sum 
vanishes for n even, as expected. We have also explicitly verified that the above expression repro¬ 
duces the value of h c * given by directly evaluating the integral in eq. (3.5) for odd values of n up to 
n = 500. Combined with the new conjecture (1.6), the above expression yields the following sum 
for the corresponding corner coefficient (for odd n > 1): 


n— 1 


cs 

0V, = 


8nn 2 (n — 1) 




k= 1 


|fc +i cos(7rfc/w) 
sin 3 (7rA;/n) 


(A.4) 


Unfortunately, the sum appearing here does not seem to be simply related to that appearing in 
eq. (3.10). 


A.2 Dirac fermion 


Now we consider the conformal dimension of the twist operator in the theory of a free Dirac fermion. 
Again the integrand in eq. (3.5) is even and so we extend the integral to the entire real line: 
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r°° du 

o 2 + sinh 2 u 

2 + sinh 2 {u/n) 

167m 2 „ 

! -00 tanh u 

lu q 

sinh u 

sinh 3 (u/n) 


(A.5) 


Again, we note that the integrand is finite at the origin, where it becomes 7(n 4 — l)/(30n) + 0{u 2 ). 
Further we note that the integrand decays exponentially fast as u —> ±oo. Next we observe that if 
we shift the contour in the complex u plane from the real axis to 3f(u) = inn, then the integrand 
is invariant up to an overall change of sign, but only for odd n. Hence we may re-express the 
conformal weight in terms of complex contour integral which runs along the contour C, illustrated 
in Fig. 7: 


1 

f du 

o 2 + sinh 2 u 
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2 + sinh 2 (w/n) 

327m 2 
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sinh u 

sinh 3 (u/n) 


(A.6) 


Hence the result can again be expressed in terms of the residues of the poles found to lie within the 
contour C. The integrand in C is only singular along the imaginary axis where sinh(u) vanishes, 
i.e., at u — ink with k = 1,2, , n — 1. Again, the first term in the integrand actually has a 

quartic and quadratic pole at each of these points and neither of these contribute to the closed 
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Table 5: Corner coefficient a n for various rational values of the Renyi index for the complex scalar 
and Dirac CFTs in d = 3. 


contour integral. Hence only the second term contributes with the result 
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Again, this result only applies for n odd. In this case, the sum does not vanish for n even, but the 
result should not be taken as the correct value of h f n . We have again explicitly verified that this 
sum reproduces the value of hJ n found by directly evaluating the integral in eq. (3.5) for odd values 
of n up to n = 500. Combining the above sum with the new conjecture (1.6), we fold the following 
expression for the corresponding corner coefficient (for odd n > 1): 


n— 1 r 


crl = 




167m 2 (n — 1) J Lsin 3 (7rA;/n) sin(7r k/n) 


(A.8) 


Here again, this new sum does not seem to be simply related to the sum in eq. (3.10) for the same 
coefficients. 


A.3 Rational values of n 

Using our integrals (3.5) for the scaling dimensions h n of the free scalar and fermion, we can 
evaluate a n in closed-form at various rational non-integer values of n. We list select results in Table 
5. Comparing the results for n = 2/3 and 3/2, we see that they explicitly satisfy the Bose-Fermi 
duality n 2 <r® c = a\j n in eq. (4.3). Similarly, comparing n = 1/6, 1/5, 1/4, 1/3 and 1/2 with the 
dual results in Table 2, we see that they obey the same duality relation. We provide a general proof 
of the duality eq. (4.3) in appendix D. 


B Sharp corner coefficient calculation 

In section 3.2.2, we computed the sharp corner coefficients K n , n = 2, 3,4 and n — oc for the complex 
scalar and the Dirac fermion. Here we explain how these calculations are performed following the 
approach of [4]. In particular, we exploit the relation between these coefficients and the so-called 
‘entropic c-functions’ for a family of corresponding massive free fields in d — 2. As was noted in 
section 3.2.2, the coefficients are given by 

<C'= where C f (mQ = t ? (B.l) 
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where S'^ s,f (m£) are the Renyi entropies of d = 2 free fields with mass m for a single interval of 
length t. The c n (t) for the complex scalar and the Dirac fermion can be evaluated as 
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& ! + <csW 


^a, cs(^ t 0) 
Ua,cs{t t CX)) 


-1 

t(log(t/2) + 2y E + (^[a] + ^[1 - a])/2) ’ 
2 

— sin(a7r)/ti_ 2a (0 , 

7T 


and 


<fW + 


w 


a, f 


W . w a ,f(t) 


+ 


l-<f (t) 




4q 2 ttq,f(t) 

t 2 !-<fW 


0, 


(B.2) 


(B.3) 


0, 

(B.4) 

(B.5) 


(B.6) 


Wa,f(£ -> 0) = —2a(log(f/2) + 2y E + (t%] + V^[-a])/2), 

2 

«a,f(£ -> oo) = - sm(an)K 2a (t). 

7T 

The differential equations can also be used to eliminate the integral in eq. (B.3) producing 

t 2 u 'l,cs - M a,cs (( 1 - 2a ) 2 + ^i 1 + <cs)) 


U?(t) = 


“lit) = 


2(1+ <J 


* 2 <f-<f(4a 2 + t 2 (l-< f )) 
2 ( 1 -< f ) 


(B.7) 


(B.8) 

(B.9) 


The final expressions for the k™' f are then given in eq. (3.22). Hence, our approach is to numerically 
solve eqs. (B.4) and (B.6) subject to the boundary conditions (B.5) and (B.7), respectively, and 
then use the results to evaluate eq. (3.22). 

Let us add here that we can evaluate K cs,i corresponding to the special case n — 1 with a similar 
approach using two-dimensional entropic c-functions. In particular, we have 

-c cs ’ f (t), (B.10) 

7r 
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with 


(Bn) 

r°° rlh 

c f (t) = 2vr/ —(B.12) 
Jo sinn [Tib) 

The values of these coefficients were evaluated by [4] and appear in Table 3 under n — 1. 


C Relation between and 


In comparing the behaviour of eq. (3.10) for large n, the authors of [18] observed a relation (3.18) 
between for the free fields and corresponding results for the logarithm of the partition function 
on the hyperbolic cylinder S 1 xi 2 , given in [28]. We will provide an explanation of this relation in 
this appendix and at the same time, we show that the same result holds for any three-dimensional 
CFT. 

Following [28], we begin by defining the logarithm of the thermal partition function evaluated 
at T = T 0 /n 

F n = — log Z (T 0 / n ), (C.l) 

which is now readily connected to the corner coefficient a n through eq. (3.11). In particular, it is 
straightforward to show that eq. (3.2) for the desired energy densities can be re-expressed as 

(G2) 

and hence eq. (3.11) can be written as 


0"n 


R 2 n 
27r Vh 2 n — 1 



(C.3) 


Now, of course, the quantity T n is infrared divergent because it involves an integration over the 
infinite volume of the hyperbolic plane HI 2 . However, both of the previous two formulae compensate 
for this divergence by dividing by Vgp to produce a finite quantity. One can produce expressions 
involving only finite quantities by focusing on the universal contribution to JF n , which amounts here 
to only retaining the finite regulator-independent contribution. The IR regulator also renders the 
volume finite and one only keeps the regulator-independent term in V jp. In particular, examining 
eq. (3.3), we have = —2nR 2 . Hence eq. (C.3) can be written as 




1 n 
47 t 2 n — 1 



d 

yJ n J n 



(C.4) 


However, before proceeding, we must first consider a difference in the conventions in the present 
paper and [28], which was alluded to in the discussion after eq. (2.7). With the heat kernel regular¬ 
ization used to produce the partition functions (3.4) for the free fields [17], one finds a nonvanishing 
energy density at T = T 0 . On the other hand, with the approach used in [28], S(T 0 ) = 0 and hence 
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d n J 7 n mv | n= i = 0 as a consequence of eq. (C.2). Since our present goal is to explain the result in 
eq. (3.18), we adopt the latter convention for the moment. Of course, this produces a simplification 
in eq. (C.4) and we are left with the following expression for the corner coefficient, 

= (&6) 

Finally, we consider the limit n —> oo for which one finds a linear growth in J r “ mv . Hence defining 
J-’oo as the coefficient of the leading term, i.e., J r “ mv ~ nT 0OJ the above expression reproduces the 
desired relation (3.18) 


a °° = 4^ jC ' 0 ° ’ (C.6) 

which was observed in [18] for the free CFTs. Of course, our analysis here made no reference to 
the free field theories discussed there and so we have established that eq. (3.18) holds generally 
for any three-dimensional CFT — assuming that the renormalization scheme is chosen such that 

S(T 0 ) = 0. 

To close, let us return to the difference in the conventions adopted here and in [28]. As noted 
above and discussed below eq. (2.7), differences in the renormalization of the Euclidean partition 
function produce a shift in T n \ 



here 




(C.7) 


Hence let us consider the form of eq. (3.18) if we instead adopt the conventions of the present 
paper. First, we note that with eq. (C.7), d n J r ™' v vanishes in the limit n —> oo. Therefore eq. (C.4) 
reduces to the following expression 


Coo 


1 

47T 2 


d F ux 

KJ n J n 


n =1 


R 3 S(T 0 ). 


(C.8) 


Again, we emphasize that although eqs. (C.6) and (C.8) are dissimilar in appearance, they both 
yield the same numerical value for a^. 


D Proof of the Bose-Fermi duality 

In this appendix, we prove the intriguing duality relation (4.3) relating the Renyi corner coefficients 
of the free scalar and the free Dirac fermion through 

n 2 < s = a\/ n , for n > 0, n G I. (D.l) 

The proof is most simply presented by adopting the conventions and notation of [28], as described 
in the previous appendix. In this case, we can use eq. (C.5) to rewrite the duality relation (D.l) as 

= 0. (D.2) 

q=l/n 


n 


Truniv, cs 
U n J ~n 


+ d q ^ 
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Now, ref. [28] provides the following simple expressions for J^mv.cs and J r “ llv,f 


3C(3) 


2 r °° 
4-/J-2 n 2 


jrnmv,cs = n ^_J_ __ dxx tanh ( A ) bg _ g—2nA) ? 


r. 


univ, f _C(3) 


= n- 


2 f 00 

^ / dX A coth(A) log (l + e~ 2nA ) 

0 


Taking derivatives with respect to n of the expressions above then yields 

9jr r ’ c ’ = |p ~ p f°<iAA 2 tanli(A) (coth(nA) - 1) . 


W"" V: * — ^A A 2 coth(A) (tanh(nA) — 1) 


(D-3) 

(D.4) 


7T 7T- 


(D.5) 

(D.6) 


We note in passing that these can be used to produce yet another integral representation of cr| s,f 
through eq. (C.5) — see eq. (7.5). Now, after some manipulations, including a change of variables 
A —y X/n in the scalar integral, we find 


Tfd n J? 




univ, f 


q=l/n 7r 


^ J n 3 -- [ dX X 2 tanh(A/n) (coth(A) — 1) , 

47T 2 TT Z J 0 

^ ^ 4 —- f dX A 2 coth(A)(tanh(A/ra) — 1). 


7T 


Hence the right-hand side of eq. (D.2) becomes 


n 3 d„F' 


cs + d t 

n 1 '-'q*' q 


univ, f 


= 03) 

3=1 /n 47T 2 


(3n 3 + 4) H —- dXX 2 [tanh(A/n) — coth A] 


7T 


(D.7) 

(D.8) 


(D-9) 


Finally, this expression vanishes identically by virtue of the exact result for the integral 
[ dX A 2 [tanh(A/n) — coth A] = — -^—^-(3n 3 + 4), V n > 0 . 


(D.10) 


This completes the proof of eq. (D.2), and hence of the duality eq. (D.l). 

If we apply the simple identity d q f(q)\ q =i/ n = —n 2 <9 n /(l/ri), we can express eq. (D.2) as 


f) x’Univ, f _ a Turov, cs 
U n' r l/n ~ uu n' r n 


(D.ll) 


which can be integrated to yield 


Jp“ v ’ f = n Jp mv,cs - / dqF f ; mv ' cs - ( JT V ’ CS - J 7 ! 


^univ, cs ^univ, f 


(D.12) 


Of course, J r , ur " v is equivalent to the universal coefficient in the sphere partition function which 
plays a role in the F-theorem [47, 56]. Eq. (D.12) shows that the thermal partition functions of the 
two free theories are not simply exchanged under n 1/n, but they are still related in a relatively 
simple manner. 
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E Relation to circular region Renyi entropies 

Let us first recall eq. (1.7): 


a„ft„L,i = 2>r (E.l) 

This relation was first observed to hold for certain holographic theories in [16] and later proven 
for general CFTs in [17]. This result was also connected to a similar relation found in [57] for 
derivatives of the Renyi entropy S n of a spherical entangling region 


d n S, 


n\ n = 1 


= —7T 2 


d +1 (d — l)T(d/2) 


T(d + 2) R d ~ l 


a 


(E.2) 


Here again, R is the radius of the sphere and V^d-i is the regulated (dimensionful) volume of the 
hyperbolic geometry HI d_1 , which appears in the construction of [23]. In particular, the following 
identity holds [17] 


r) <? I 
u n k -'n\ n= \ 


d — 1 Vjjd—i 

2 R d ~ l 


f) h I 

u n' L n\ n= i 


(E.3) 


and more generally, 14 for k > 1 



n= 1 


d — 1 Vjjd-i 

k + 1 R d ~ l 


yi (-) m k\ 
4^ (k — m)\ 

m— n V 7 


pjk—m 

'-'n 


hr, 


I n=1 


(E.4) 


Note that the above sum does not extend to m = k because hi = 0. The relations between these 
derivatives can be written in a more compact form as 



n =1 


R 


d -1 


{d — 1) V B d-i 


(k + i)d k n s n \ n=1 + e d k n -%\ n=1 


(E.5) 


where we drop the second term for k = 1. 

Note that the scaling dimension on the left-hand side is a finite quantity of this last identity 
(E.5), while the Renyi entropies on the right-hand side are UV divergent. Of course, this divergence 
is compensated by dividing by I4 p to produce a finite quantity. One can produce an identity relating 
finite quantities by focusing on the universal contribution to the Renyi entropy, which is either the 
coefficient of the logarithmically divergent term in an even number of dimensions or the finite 
contribution in an odd number of dimensions. For example, let us focus on d = 3 where we can use 
eq. (3.3) to extract the finite term from Vh 2 - Then eq. (E.5) yields 



71—1 


1 

47T 


(k + 1 ) d k n Sl 


2 ofc—1 rmmv I 


171=1 


+ k d k ~ L S, 


I 71=1 


(E.6) 


It is straightforward to explicitly verify this relation in, e.g., the holographic model in section 5. 
14 Please note that the original expression for k > 2 in eq. (2.53) of [17] was incomplete. 
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Now, expanding our new conjecture (1.6) around n — 1 yields 


d n h n \ n= i = tt k d n a n \ n =i , 
for k > 1. Hence we can combine this result with eq. (E.6) as 


d k a I = — 

" "ln=l 4^2 


fc + 2 
fc + 1 


ar'sTL, +(* + 1 ) D's, 


k rmniv I 


™ ln=l 


(E.7) 


(E-8) 


which is valid for k > 1. However, eqs. (E.6) and (E.7) yield an additional relation for k — 0, which 
can be seen as another realization of our original conjecture (1.5): 


Qumv 

°n D n 


7=1= 27r2(J l- 


(E.9) 
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